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Abstract 

On the basis of features observed in the exact perturbation approach solution for the 
eigenspectrum of the dilute A3 model, we propose expressions for excitations in the 
dilute A^ and Aq models. Principally, we require that these expressions satisfy the 
appropriate inversion relations. We demonstrate that they give the expected Ej and 
Eq mass spectra, and universal amplitudes, and agree with numerical expressions 
for the eigenvalues. 



1 Introduction 



The dilute Al model is an exactly solvable, restricted solid-on-solid model 
denned on the square lattice. At criticality, the model can be constructed [1,2] 
from the dilute 0(n) loop model [3,4]. Each site of the lattice can take one of 
L possible (height) values, subject to the restriction that neighbouring sites 
of the lattice either have the same height, or differ by ±1. Most importantly, 
the model can also be solved away from criticality. The off-critical Boltzmann 
weights of the allowed height configurations of an elementary face of the lattice 
are parametrised in terms of elliptic theta functions [1] . The interpretation of 
the elliptic nome p differs according to whether L is even or odd. In particular, 
for L odd the up-down symmetry of the Boltzmann weights is broken away 
from criticality. For L = 3 the elliptic nome plays the role of magnetic field. 
Moreover, the dilute A 3 model provides, in one of its regimes, an integrable 
lattice realisation of the Eg Ising model, being in the same universality class 
as the two-dimensional Ising model in a magnetic field. 

The calculation of the various off-critical thermodynamic properties of the 
model have verified this correspondence. The singular part of the bulk free 
energy of the dilute A 3 model in the appropriate regime gives the magnetic 
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Ising exponent 5 = 15 [1], which also follows from the calculation of the local 
height probability [5] . The expected Ising magnetic surface exponent 5 S = — y 
follows from the excess surface free energy [6] . Moreover the E 8 mass spectrum, 
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predicted by Zamolodchikov [7,8] for the Ising model in a magnetic field is 
seen in the single particle excitation spectrum [9-12]. Here the masses are 
normalized such that m\ — 1. They coincide with the components of the 
Perron- Frobenius vector of the Cartan matrix of the Lie algebra Eg. 

In this paper we consider off-critical excitations in the dilute A 4 and A 6 models, 
which are expected to be related to the E 7 and E 6 scattering theories. The E 6 
masses are (see, e.g., [13-15] and refs therein) 

mi = mi — 1 

m 2 = 2cosf =1.414 213... 

m 3 = mi = 2 cos ^ = 1.931 851 .. . 
m 4 = 4cosfcos^ = 2.732 050... 

The E? masses, with m± = 1, are [13-15] 

m 2 = 2cosf| = 1.285 575... 

m 3 = 2cosf = 1.879 385... 

m 4 = 2cosf| = 1.969 615... 

m 5 = 4 cos f| cos ff = 2.532 088 .. . 
m 6 = 4 cos f cos ^ = 2.879 385 .. . 
m 7 = 4 cos f| cos f = 3.701 666... 



(1.3) 



Our approach begins in the next two sections by considering the inversion 
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relations that hold for the off-critical dilute Al models, how our solution [12] 
satisfies them in the case L = 3, and how the E s structure manifests itself 
within the solution. In the subsequent sections we propose solutions for A A 
and A 6 and demonstrate the expected E 7 and E 6 mass spectra. We conclude 
with some numerical evidence and discussion. 



2 Inversion relations 



The eigenvalues of the row transfer matrix of the dilute Al model, defined on 
a periodic strip of width N, where we take N even, are [9] 
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where the N roots u,- are given by the Bethe equations 
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(2-1) 
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with a; = exp(i7r£/(L + 1)) for £ = 1,...,L. For regime 2, the regime to 
be considered, the spectral parameter u lies in the range < u < 3A, with 
A = ns/r, where s = L + 2 and r = 4(L + 1). 

The standard elliptic theta functions 0i(m), 04 (u) of nome p are defined as 



0i( M )=0 1 ( M ,p) = 2p 1 / 4 shiM J] (l-2p 2n cos2u + p 4n ) (l-p 2n ) , (2.3) 

n=l 

oo 

4 (u) = 4 (u, p) = II ( x - 2 ^ 2n_1 cos 2u + P An ~ 2 ) i 1 ~ P 2n ) ■ ( 2 - 4 ) 



n=l 



Also of use are the conjugate variables 

w = e -2*u/e a . = e -r»/« j 



(2.5) 
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where no me p = e e . The relevant conjugate modulus transformations are 
& 1 (u,p) = (^y 2 e-<- u -* /2) * /e E(w,q 2 ), (2.6) 

/ 7r \l/2 

4 («,P) = (^) e-^Wf'Ei-wrf), (2.7) 
where g = e _7r /<E and 

oo 

= n a - - p n ^ i ){\ - p n ). (2.8) 

71=1 

For this model, the partition function per site k was first calculated using the 
inversion relation [1,5] 

k(u) k(u + 3A) - ^ l(2A = ^ l(3A = ^ l(2A + ^ l(3A + u) (2 9) 
k(u) k[u + S\) - <5>?(2A)0?(3A) • { } 

In this way the bulk free energy per site / = log k was found to be 

k{\ - x 2 - k ){l + x^ sk ) ' (2 ' 10) 

The same result was derived [12,16] from the Bethe Ansatz solution for the 
groundstate eigenvalue A (u). 

Making use of the Poisson summation formula in the free energy (2.10) the 
leading singularity as p — > in regime 2 is 

f~Ap r ^ s , (2.11) 

where the amplitude A is given in terms of L by 

cos W+en 

A- U /o C0S Ul +2 )J (212) 

sin r 2 " (L+i) v 1 j 

V 3(L+2) ) 



and we have taken the isotropic value u = 3A/2. 

Excitations in the eigenspectrum can be considered in terms of the quantity 

rj{u) = lim -44. 2.13 

Jy ' N-+oc An(u) V ' 
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The inversion relation (2.9) is simply 

Tj{u) r-j{u + 3 A) = 1, 



(2.14) 



but there is a further relation to be satisfied [12], 

r j (u)r j (u + 2\)=r j (u + \). (2.15) 

Our approach here is not to solve the inversion relations directly, as was done, 
e.g., by Kliimper and Zittartz for the excitation spectra of the eight- vertex 
model [17]. Rather, in the light of our results for the excitations of the dilute 
A3 model, we use the above inversion relations to give constraints on the Lie 
algebraic properties of a conjectured solution. We then test our results as best 
we can by numerically diagonalising the transfer matrix, and by comparison 
with results for E 7 and E e obtained by other methods. 



3 The dilute A% model and the E% mass spectrum 



We now summarise our results [12] for the dilute A3 model, obtained by the 
exact perturbation approach [18]. The leading excitations in a given band of 
eigenvalues can be written in the compact form 

r (w) = TT Et^hl ^°)E(-^- a /w, x™) 
A > l } E{-x a w, x™)E{-x™- a w, x 60 ) ' 1 ; 



where the numbers a and n(a) are given in Table 1. The E 8 numbers a have 
been discussed by McCoy and Orrick for the related Hamiltonian [11]. They 
appear, e.g., in E 8 scattering theory [14] and in E 8 Lie algebraic polynomials 
[19]. The number n(a) denotes the relevant band of eigenvalues. 



Note that within a band of eigenvalues there may be more than one class 
of excitation. For example, in the leading band of eigenvalues there are two, 
which arise from a 2-string and a 4-string structure in the Bethe roots [9,10]. 
The expression (3.1) is the leading excitation for each class of eigenvalue. The 
last excitation within a class is also given by (3.1), but with positive argument 
in the elliptic functions. 

In the original variables (3.1) reads 



rAu 
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Table 1 

Parameters appearing in the eigenvalue expression (3.1) 
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4 
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8 


6 
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The various correlation lengths follow as 

ZJ 1 = -logrj(u), 



(3.3) 



where we take the relevant leading eigenvalue at the isotropic point u = 3 A/2, 
which for L = 3 is u = 

The fundamental correlation lengths can thus be written 

Ufo + hP 8/15 )M^ + hP 8/15 ) 
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Expanding this expression in powers of p gives 

rrij ~ 8p 8/15 ^sinff as p^O. 



(3.5) 



This is the formula obtained by McCoy and Orrick [11] for the Hamiltonian, 
from which the E$ masses in (1.1) are recovered by virtue of trig identities. 



In particular, 



ff/ 1 ~8p 8/15 (sinfr + sin^f) = 16 sin f cos § p 8/15 



(3.6) 
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We are now able to consider the universal magnetic Ising amplitude [16]. From 
(2.11) and (2.12), 



/^^^p 16 / 15 as p^O. (3.7) 
cos ^ 



Combining this with (3.6) gives 

fil = =— ^ = 0.061 728 589... as p^O. (3.8) 

16 v 3 sin | cos ^ 



This is the result for the universal magnetic Ising amplitude obtained earlier by 
thermodynamic Bethe Ansatz calculations based on the E 8 scattering theory 
[15] (see also Ref. [20] in the context of the form-factor bootstrap approach). 
Here it has been obtained from the lattice model. 

From the outset, no assumptions were made on the E 8 structure in the dilute 
A 3 model, both in our own calculations, and in the thermodynamic Bethe 
Ansatz calculations [9,11]. We now highlight a few of the E 8 features as a 
guide to our considerations of E 7 and Eq. 

First, each a value occurs in (3.1) together with its complement in 30, the 
Coxeter number of E s , so that no integer greater than 15 appears in the lists 
in Table 1. 

Second, the inversion relation 

r j {w)r j {x 30 w) = l, (3.9) 



is satisfied by an expression of the form (3.1) for any a value. However, the 
stronger relation 

rj (w) rj (x 20 w) = rj (x w w) , (3.10) 



is satisfied if, within the set of integers, one finds not only a, where a = 1, . . . , 9, 
but also a + 10, or equivalently its complement in 30, 20 — a, by virtue of the 
properties 

E(z,p) = E(p/z,p) = -zE(z-\p). (3.11) 



The integer a = 10 may appear alone, because the factor it contributes to 
r-j{w) satisfies (3.10) by itself. From Table 1, the sets of integers found by the 
perturbative approach [12] all have these features. 
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Finally, we observe that the nome p cancels in (3.8) because of the relationship 
between the power of p occurring in / and in £ tl . Indeed, this combination 
defines the hyperscaling relation between the corresponding critical exponents. 



4 The dilute A 4 model and the E 7 mass spectrum 



We now use our observations for E$ to arrive at a conjecture (equivalent to 
(3.1)) for the excitations of E 7 . 



The free energy expression is, from (2.11) and (2.12), 

f-^lp 10/9 ^ p^O. (4.1) 



sin lg 



In order to obtain a finite expression from f £f, we thus expect 

^r 1 ~ p 5/9 as p -> . (4.2) 

This power of the nome must appear in the expression equivalent to (3.2) for 
E 7 , and is thus related to the one we must propose for Tj(w) by the conjugate 
modulus transformation (2.7), namely 

where for L = 4, x = e _7r ^ 20e . 

The inversion relation in conjugate modulus form is 

r j (w)r j (x 36 w) = 1. (4.4) 



However, the Coxeter number for E 7 is 18, so that we expect to select our 
integers from 1, . . . , 9. We thus propose that the excitations for E 7 take the 
form 

r (w) - n Et^hh * 72 )E(-x^/w, x^) 
]K 1 \ L E{-x 2a w, x 72 )E{-x 36 - 2a w, x 72 ) ' V ' } 



The additional relation which serves to constrain the possible a values is 

r j( w ) rj(x 24 w) = r-j(x 12 w). (4.6) 
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Table 2 

Parameters appearing in the eigenvalue expression (4.5). 
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3, 5, 7, 9 



This condition is satisfied if, within a set of possible integers, a appears to- 
gether with a + 6 or equivalently 12 — a, apart from a = 6 whose contribution 
satisfies (4.6) by itself. 



Written in terms of the original variables the expression (4.5) is 

r(u) TJ ^-^P^M^-f,^ 9 ) (A7) 



Taking the isotropic value u = 9n/20 we obtain 

m r' 2ri° / 4(S + f ' p5/8> 



for the masses, and so 

rrij ~ 8 p 5/9 sin if as P • ( 4 - 9 ) 

a 



We now turn to the sets of integers associated with E-j in the context of Lie 
algebraic polynomials [19] which form the first six rows of the third column 
of Table 2. Clearly these integers satisfy the constraints described above as 
being placed upon them by (4.6). Together with the last row, they are also to 
be found within the table given for E 7 scattering in [14] . 



Applying trig identities to the sum in (4.9) with these sets of integers gives 
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which, apart from normalisation, correspond to m\, . . . , of (1.3). Q As an- 
other piece of evidence for our identification of a = 6 with m\, from which the 
others follow, we consider the amplitude 



2 2^ 1 _ 1 

sing ' (8sinf) 2 ~ 8V3cos~ T 



Id -^-7——^ — ^ 57 ^ p^O. (4.11) 



This is in agreement with the E 7 thermodynamic Bethe Ansatz result [15]. 



5 The dilute A§ model and the E 6 mass spectrum 



The free energy expression for the dilute A§ model is, again from (2.11) and 
(2.12), 

f-^TP 7/6 as p-+0, (5.12) 
cos ■ 



12 



and so we expect 



^~// 12 as (5.13) 



1 There is another relationship between the E-j mass ratios, the trigonometric 
expression of (4.9) and integers in the table of [14]. However, the one described here 
is necessary in the context of the solvable dilute model in order to satisfy its 
inversion relations. 
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This power of the nome must appear in the expression equivalent to (3.2) 
for Eq, and thus prescribes the nome of the expression we propose for Tj(w), 
because in the conjugate modulus transformation (2.7), 

e -*/i2 - e" 24 ^ = x 96 , (5.14) 

where in the case L = 6, x = e"" 2 / 28e . The inversion relation in conjugate 
modulus form is 

r j( w ) rj(x 48 w) = 1. (5.15) 



Finally, the Coxeter number for E Q is 12, so that we expect to select our 
integers from 1, . . . , 6. We thus propose that the excitations for E 6 take the 
form 

A >~ \ l E(-x^w,x^)E(-x^w,x^) ■ 



The additional relation which serves to constrain the possible a values is 

r j (w)r j (x 32 w) =rj(x 16 w). (5.17) 

Thus within any set of possible integers, a must appear together with a + 4 or 
equivalently 8 — a, apart from a = 4 which satisfies (5.17) by itself. Written 
in terms of the original variables the expression (5.16) is 

A jaw _ 7u V 7/12\M / (12-a)7T _ Ju 7/12N 
r ( M ) = ]- r 4 ^4 12»P ^24 12»P ^ (5 . 18) 

JW a ^(l + f|,P 7 / 12 K4( ii ^ + fi ! p 7 / 12 ) 

Taking the isotropic value u = 3ir/7 we obtain 

,-r=2E < (5.19) 



m 

24 4 ' 



for the masses. Thus 

rrij ~ 8 p 7/12 2 sin ff as p -> . (5.20) 



The integers given in Table 3 satisfy the constraint placed upon them by 
(5.17). Apart from the entry for j = 4, these integers are again to be found in 
[19], and they appear within the table of [14] for E 6 . 
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Table 3 

Parameters appearing in the eigenvalue expression (5.16). 
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n{a) 


a 


1 ,1 


1 


4 


2 


2 


1, 5 


3, 3 


2 


3, 5 


4 


3 


2, 4,6 



Applying trig identities to the sum in (5.20) with these sets of integers gives 

£ sin ff = V3/2, 
^ sin f = v^/v^, 

o=l,5 

X] sin ff = v / 3cos^, (5.21) 

a=3,5 

^ sin ff = ^008^, 

a=2,4,6 

which, apart from normalisation, correspond to mi, ... , m 4 of (1.2). Our iden- 
tification of a = 4 with mi, gives the amplitude 

f^l = ^..— }— = ——}: as p ^o, (5.22) 

^ COS iL ( 4v ^)2 2^(1 + ^) 

which is in agreement with the thermodynamic Bethe Ansatz result [15]. 



6 Numerical evidence and discussion 

Based on our result (3.1) for the eigenspectrum of the dilute A 3 lattice model in 
regime 2, and its resulting E 8 structure, we have proposed analogous formulae 
for the dilute A 4 and A$ models under the assumption of corresponding E 7 and 
E 6 structures. Such correspondence is to be expected on a number of grounds. 
For example, at criticality the central charges of the dilute A L models are 
known from the underlying loop model [1]. In regime 2, c = ^ for the A A 
model and c = | for the A 6 model. These are the same as the E 7 and E 6 
values [13]. 

A number of considerations have motivated our final results. Our first input 
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was the hyperscaling relation, /£ 2 = constant, which constrains the power of 
the elliptic nome p appearing in the inverse correlation lengths. We found that 
the stronger inversion relation (2.15) constrains the set of integers a appearing 
in the eigenvalue expressions. We took these numbers from the Lie algebraic 
polynomials tabulated by Kostant [19]. Our results produce the E e (1.2) and 
Ej (1.3) masses in the critical limit p — * 0. However, the configuration of a's 
for the heaviest mass does not appear in the Kostant polynomials. We chose 
that configuration to be consistent with the predicted E e and E 7 mass spectra, 
and subsequently noted that it had been observed in the context of scattering 
theory [14]. Our identification of the a's associated with the lightest masses 
also gives the universal amplitudes (5.22) and (4.11) in agreement with the 
thermodynamic Bethe Ansatz results based on the E 6 and E 7 algebras [15]. 

We have performed a number of numerical tests on the eigenspectra of the 
dilute A A and A Q models to check our results. Specifically, we have diagonalised 
the periodic row-transfer matrix for finite lattice sizes. Consider the dilute 
Aq model first. Here A = ^F. The largest eigenvalue Ao is 3-fold degenerate 
in the thermodynamic limit. Successive numerical estimates of the first few 
eigenvalue ratios A /Aj at the isotropic point u = 3 A/2 are tabulated in Table 
4 for the values p — 0.1 and p = 0.3. Excellent agreement is seen with the 
expected result (5.18), which reduces to 



is.-, n 



(6.1) 



The dilute A 4 model at A = ^ is more complicated. A detailed numerical 
study of the Bethe Ansatz equations has revealed all seven masses [21]. How- 
ever, the eigenvalue spectrum is dependent on the sign of p. In this case, all of 
the E-j masses are observed in the p < regime (regime 2~). Only a subset is 
observed for p > (regime 2 + ). Our numerical results for the first few leading 
eigenvalues are shown in Table 5 for p = —0.3. The eigenvalues Ai and A 3 
do not appear in the eigenspectrum for p = 0.3. Clearly there is excellent 
agreement with our result (4.7), which here simplifies to 



Ao 
A,- 



n 



04 



36 

- 36 



f,p 5/9 ) 

f,pV9) 



(6.2) 



We expect this result to hold in regime 2~ for all of the masses, or correspond- 
ingly for each set of a's given in Table 2. Apart from A 1 and A 3 , we have not 
explored further which of the eigenvalues are absent in regime 2 + . We await 
the publication of Ref. [21]. 

In contrast with the dilute A A model, the mass spectrum of the dilute A 6 model 
appears to be equivalent in regimes 1 ± . Such equivalence holds for the dilute 
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Table 4 

Numerical estimates with increasing system size N of leading eigenvalue ratios in 
the dilute A§ model at A = Also shown is the expected exact result (6.1) in the 
thermodynamic limit. The corresponding values of a are as given in Table 3. 





N 


Ao/Ai 


A /A 2 


A0/A3 


p = 0.1 


3 


6.0279 








4 


6.7882 








5 


6.9281 








6 


6.9474 


15.268 






7 


6.9501 


15.511 


41.05 




oo 


6.9505 


15.590 


42.44 


p = 0.3 


3 


89.93047 








4 


90.08438 








5 


90.08605 


652.6278 






6 


90.08607 


652.7399 


6434.75 




7 


90.08607 


652.7442 


6436.87 




oo 


90.08607 


652.7444 


6437.08 



A L models with L odd, where the eigenspectrum is independent of the sign of 
p. This is a consequence of the off-critical weights breaking the Z 2 symmetry 
for L odd. However, for L even this symmetry is not broken. As to why the 
mass spectrum may be the same for the dilute A e model in regimes 2 ± , this 
remains one of the mysteries of the dilute Ai models for L even, which are 
yet to be investigated. 

Finally we note that although the evidence for our conjectured results is con- 
vincing, they of course await a formal derivation. 
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Table 5 

Numerical estimates with increasing system size N of leading eigenvalue ratios in 
the dilute A<± model at A = jq. Also shown is the expected exact result (6.2) in the 
thermodynamic limit. The corresponding values of a are as given in Table 2. 





N 


Ao/Ai 


A /A 2 


A0/A3 


A0/A4 


p = -0.3 


4 


116.09490 


492.5475 








5 


116.09969 


493.2263 


8669.13 






6 


116.09973 


493.2413 


8724.17 


11928 




7 


116.09973 


493.2416 


8726.53 


12067 




oo 


116.09973 


493.2416 


8726.64 


12190 
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